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Abstract — We consider the problem of controlling the prop- 
agation of an epidemic outbreak in an arbitrary contact 
network by distributing vaccination resources throughout the 
network. We analyze a networked version of the Susceptible- 
Infected-Susceptible (SIS) epidemic model when individuals in 
the network present different levels of susceptibility to the 
epidemic. In this context, controlling the spread of an epidemic 
outbreak can be written as a spectral condition involving the 
eigenvalues of a matrix that depends on the network structure 
and the parameters of the model. We study the problem of 
finding the optimal distribution of vaccines throughout the 
network to control the spread of an epidemic outbreak. We 
propose a convex framework to find cost-optimal distribution 
of vaccination resources when different levels of vaccination 
are allowed. We also propose a greedy approach with quality 
guarantees for the case of all-or-nothing vaccination. We 
illustrate our approaches with numerical simulations in a real 
social network. 

I. Introduction 

Motivated by the problem of epidemic spread in human 
networks, we analyze the problem of controlling the spread 
of a disease by distributing vaccines throughout the indi- 
viduals in a contact network. The problem of controlling 
spreading processes in networks appear in many different 
settings, such as epidemiology JTJ, 0, computer viruses 
0, or viral marketing [4J. The dynamic of the spread 
depends on both the structure of the contact network, the 
epidemic model and the values of the parameters associated 
to each individual. We model the spread using a recently 
proposed variant of the popular SIS epidemic model in 
which the infection rate is allowed to vary among the set 
of individuals in the network [5|. In our setting, we can 
modify the individual infection rates, within a feasible range, 
by injecting different levels of vaccination in each node. 
Injecting a particular level of vaccination in a node has 
also an associated cost, which can vary from individual 
to individual. In this context, we propose efficient convex 
framework to find the optimal distribution of vaccination 
resources throughout the networks. 

The dynamic behavior of spreading processes in networks 
have been widely studied. In J6[, Newman studied the epi- 
demic thresholds on several random graphs models. Pastor- 
Satorras and Vespignani studied viral propagation in power- 
law networks 0. This initial work was followed by a long 
list of papers aiming to study the spread in more realistic 
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network models. Boguna and Pastor-Satorras |8| considered 
the spread of a virus in correlated networks, where the 
connectivity of a node is related to the connectivity of its 
neighbors. In (9), the authors analyze spreading processes 
in random geometric networks. The analysis of spreading 
processes in arbitrary contact networks was first studied by 
Wang et al. ifTOl for the case of discrete-time dynamics. 
In ifTTI . Ganesh et al. proposed a continuous -time Markov 
process to relate the speed of spreading with the largest 
eigenvalue of the adjacency matrix of the contact network. 
The connection between the speed of spreading and the spec- 
tral radius of the network was also found for a wide range 
of spreading models in fl2l . The relationship between the 
spectral radius of a contact network and its local structural 
properties were explored in [13], |14|. 

The development of strategies to control the dynamic of 
a spread process is a central problem in public health and 
network security. In JT5 |, Borgs et al. proposed a proba- 
bilistic analysis, based on the theory of contact processes, 
to characterize the optimal distribution of a fixed amount 
of antidote in a given contact network. In [16], Aditya et 
al. proposed several heuristics to immunize individuals in a 
network to control virus spreading processes. In the control 
systems literature, Wan et al. proposed in [17| a method 
to design optimal strategies to control the spread of a virus 
using eigenvalue sensitivity analysis ideas together with con- 
strained optimization methods. Our work is closely related 
to the work in [18| and [19|, in which a continuous-time 
time Markov processes, called the N-intertwined model, is 
used to analyze and control the spread of a SIS epidemic 
model. 

In this paper, we propose a convex optimization frame- 
work to efficiently find the cost-optimal distribution of 
vaccination resources in an arbitrary contact network. In our 
work, we use a heterogeneous version of the N-intertwined 
SIS model Q to model a spread process in a network of 
individuals with different rate of being infected and recov- 
ered. We assume that we can modify the rates of infection of 
individuals, within a feasible range, by distributing vaccines 
to the individuals in the network. We assume that there is a 
cost associated to injecting a particular amount of vaccina- 
tion resources to a each individual, where the cost function 
can vary from individual to individual. Our aim is to find the 
optimal distribution of vaccination resources throughout the 
network in order to control the spread of an initial infection 
at a minimal cost. We consider two version of this problem: 
(;) The fractional case, in which we are allowed to inject 



a fractional amount of vaccination resources in each node 
of the network, and (ii) the combinatorial case, in which 
we fully vaccinate a selection of individuals in the network, 
leaving the rest of nodes unvaccinated. 

The paper is organized as follows. In Section II, we 
introduce our notation, as well as some background needed 
in our derivations. In Section III, we formulate our problem 
and provide an efficient solution based on convex optimiza- 
tion. In Section IV, we study a combinatorial version of 
the problem studied in Section III and provide a greedy 
heuristic algorithm with a quality guarantee. We include 
some conclusions in Section V. 

II. Notation & Preliminaries 

In this section we introduce some graph-theoretical 
nomenclature and the dynamic spreading model under con- 
sideration. 



for i — 1, . . . , n. 

As proved in Q, the exact probability of infection is 
upper bounded by its approximation pi (t). A local stability 
analysis of the above system of ODE's around the disease- 
free equilibrium, pi = for all i, provides the following 
result 0: 

Proposition 1: Consider the N-intertwined SIS epidemic 
model in (QJ. Then, an initial infection converge to zero 
exponentially fast if 

Ai (A g )< 6 -. 

The above provides a simple condition to guarantee a con- 
trolled epidemic dynamics in terms of the largest eigenvalue 
of the adjacency matrix. In the following, we derive a 
similar condition when the infection parameters vary from 
individual to individual within the network. 



A. Graph Theory 

Let Q — (V,£) denote an undirected graph with n 
nodes, m edges, and no self-loop^] We denote by V (G) = 
{vx,..., v n } the set of nodes and by £ [Q) C V (G) x V (G) 
the set of undirected edges of G- If {i,j} £ £ (G) we call 
nodes i and j adjacent (or neighbors), which we denote by 
i ~ j. We define the set of neighbors of a node i e V 
as N t = {j £ V{G) : {i,j} £ £(G)}- The number of 
neighbors of i is called the degree of node i, denoted by di. 
The adjacency matrix of an undirected graph G, denoted by 
Ag = [o>ij], is an n x n symmetric matrix defined entry- 
wise as Ojj = 1 if nodes i and j are adjacent, and = 
otherwise^] Since Ag is symmetric, all its eigenvalues, 
denoted by \\{Ag) > \2(Ag) > . . . > X n (Ag), are real. 

B. N -Intertwined SIS Epidemic Model 

Our modeling approach is based on the N-intertwined SIS 
model proposed by Van Mieghem et at. in 0. In contrast 
with previously proposed models, the N-intertwined model 
is a continuous-time networked Markov process with 2™ 
states able to model the dynamics of a viral infection in an 
arbitrary contact network. Using the Kolmogorov forward 
equations and a mean-field approach, one can approximate 
the dynamics of the viral spread using a system of n ordinary 
differential equations, as follows. Consider a network of n 
individuals described by the adjacency matrix Ag = [a^]. 
The infection probability of an individual at node i € V (G) 
at time t > is denoted by Pi(t). Let us assume, for now, 
that the viral spreading is characterized by two positive 
parameters-a constant infection rate j3 > and a curing 
rate 6 > 0. Hence, the N-intertwined SIS model in J3) is 
described by the following set of n ODE's: 



dpi jt) 
dt 



= (1- Pi (*)) p a vPi W - 5 Pi (*) ' C 1 ) 



An undirected graph with no self-loops is also called a simple graph. 
2 For simple graphs, an = for all i. 



C. Non-Homogeneous N-Intertwined SIS Epidemic Model 

A direct extension of the N-intertwined model for node- 
specific infection and curing rates, and 5i, is 

dpi (t) 



dt 



= (1 - pi (t)) ^ ^ aijpj (t) - SiPi (t) . 

3=1 



We can write the above dynamics in matrix form as 
dp(t) 



dt 



(BA g -D)p(t)-P(t)BAgp(t), (2) 



where p(t) = (p x (t) , . . . ,p n (i)) , B = diag(Pi), D = 
diag(Si), and P (t) — diag(pi). Concerning the non- 
homogeneous epidemic model, we have the following result: 

Proposition 2: Consider the heterogeneous N-intertwined 
SIS epidemic model in Q. Then, if 

Ai (BA — D) < -e, 

an initial infection p (0) £ [0,1]™ will converge to zero 
exponentially fast, i.e., there exists an a > such that 

\\Pi (t)\\ < a \\pi (0)|| e' et , for all t > 0. 
Proof: First, we have 



dpi (t) 
dt 



= j3i ^ aijP 3 (t) - SiPi (i) - PiPi (t) 2J aijPj (*) 



3=1 



3=1 



<Pi ^ a ^Vj (t) - SiPi (*) , 

3=1 

since pi, Si, pi (t),a,ij > 0. Therefore, the linear dynamic 
system 



dpi (t) 
dt 



Pi Y a vP3 ( f ) _ S iPi (*) ' 



(3) 



3=1 



upper-bounds the nonlinear dynamical system Q when they 
share the same initial conditions, i.e., p (t) > p (t) for t > 
when p (0) — p (0). 



This linear dynamic system can be written in matrix form 



x 10- 3 



as 



dp(t) 
dt 



= (BA g -D)p(t). 



For the above linear system to be stable, we need the 
eigenvalues of BA — D to be in the open left half-plane. The 
state matrix BAg — D has real eigenvalues, since it can be 
transform via a similarity transformation to the symmetric 
matrix B 1 / 2 AgB 1 / 2 — D. Hence, exponential asymptotic 
stability, with an exponential rate e, is equivalent to the 
largest eigenvalue Ai {BAg — D) < —e. ■ 
In the above analysis, we have shown that the linear dy- 
namics in <|3j upper-bounds the mean-field approximation in 
(|2]); thus, the spectral result in Proposition [2] is a sufficient 
condition to control the evolution of an epidemic outbreak. 
In the following section, we use this result to characterize 
the profiles of infection rates that results in a stable linear 
dynamics. 

III. A Convex Framework for Optimal Resource 
Allocation 

Our main aim is to propose an efficient optimization 
framework to find the optimal distribution of vaccines to 
control the spread of an epidemic outbreak in a given net- 
work. In this section, we consider the fractional vaccination 
problem. In the fractional case, we assume that we are 
able to modify the infections rates ft in the network by 
distributing vaccination resources throughout the individuals 
in the network. We assume that the infection rates of each 
individual can be modified within a particular feasible inter- 
val, ft. < ft < Pi, where ft > is the value of the natural 
infection rate for node i, which is achieved in the absence 
of any nodal immunization, and ft > is the minimum 
possible infection rate for node i, which is achieved when 
we allocate a large amount of vaccines at node i. In Section 
|IV[ we will consider a combinatorial version of the above 
fractional strategies. In the combinatorial case, we will 
assume that the infection rate can only take one of two 
values, Pi <E ^P^Pi^. In the fractional case considered in 
this section, we propose an optimization framework to find 
the optimal distribution of resources when there is a cost 
function function associated to different values of Pi. 

A. Vaccination Cost 

The cost of achieving a particular infection rate for node 
i is denoted by fi (Pi). This cost function is node-dependent 
and presents the following properties: 

1) The cost of achieving the natural infection rate is zero, 
i.e., fi (ft) = 0. 

2) The maximum cost of vaccinating node i, denoted by 
Tj, is achieved at the minimum infection rate, i.e., 
max^j. fi (Pi) = fi (pA = Ti. 

3) The vaccination cost function is monotonically de- 
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Fig. 1. Convex cost function in J3J. 



Apart from the above properties, we make the following 
convexity assumptions on the cost function /,; to obtain a 
tractable convex framework: 

Assumption 1: The vaccination cost function, fi (Pi), is 
twice differentiable and satisfies the following constrain: 



ft {Pi)>-jfi(Pi), 



(4) 



for pi e [ft., ft 

Notice that, since fi is monotonically decreasing, we 
have that f[ (Pi) < 0; thus, we have that Assumption 1 
implies that /" (Pi) > 0. In other words, Assumption 1 is 
stronger than convexity. For example, a function that satisfies 
Assumption 1 with equality is: 



fi (Pi) = % 



fir 1 -ft 



(5) 



creasing in the interval Pi g 



In practice, for low values of P i and Pi, this function takes a 
shape of practical interest. For example, in Fig. [T] we plot the 
function in ^ for §_. = 1.75e-3, ft = 8.66e-3, and T t = 
1. In the abscissa of this plot, we represent the vaccination 
cost fi (Pi), which is in the range [0, 1]. We observe how 
the cost function is convex and presents diminishing returns, 
since the reduction in the infection rate for a given amount 
of investment is greater in the low-cost range than in the 
high-cost range. 

B. Problem Statements 

In this subsection we propose an optimization framework 
to find the cost-optimal allocation of vaccines in a given 
contact network Q with adjacency matrix Ag. In particular, 
we consider the following problem: 

Problem 1: Given a curing rate profile, {Si : i G V (Q)\, 
and a vaccination cost function fi (Pi) for Pi € P^Pi , 
find the optimal distribution of vaccines to control the 
propagation of an epidemic outbreak with an asymptotic 
exponential decaying rate e at a total minimum cost. 



According to Proposition [2] this problem can be mathemat- 
ically stated as the following optimization problem: 

T* = min {fs%} Yn=x ft (A) 

s.t. Ai (BAg -£>)<-£ (6) 

i = l,...,n, 

In the following subsection, we propose a convex formu- 
lation to solve this problem under Assumption 1. 

C. Semidefinite Programming (SDP) Approach 

Our formulation is based on writing the spectral stability 
condition Ai (BAg — D) < — e using a simple semidefinite 
constrain. In particular, we have the following result: 

Lemma 3.1: For Ag symmetric, B — diag (Pi), and D = 
diag (Si), we have that Ai (BAg — D) < — e if and only if 
(D - el) B- 1 -Ag>Q. 

Proof: Notice that BAg — D is a matrix similar to 
B 1 / 2 AgB 1 / 2 — D, since we can pre- and post- multi- 
ply the former matrix by B~ x / 2 and B 1 ! 2 , respectively, 
to obtain the latter. Hence, since B 1 ! 2 AgB 1 / 2 - D is 
a symmetric matrix with real eigenvalues, the eigenval- 
ues of BAg — D, including Ai (BAg), are all real. 
Then, we have that A x (BAg — D) < —e if and only if 
A, ((D - el) - BAg) = A, ((£> - si) - B^AgB 1 / 2 ) > 
0, which is equivalent to (D - el) - B 1/2 A g B 1 / 2 > 
0. Applying a congruence transformation to (D — el) — 
B 1 / 2 AgB 1 / 2 by pre- and post-multiplying by B~ x / 2 , 
we obtain that Ai (BAg — D) < — e if and only if 
(D - el) B- 1 - Ag y 0. ■ 
Using the above Lemma, we can rewrite the optimization 
problem [T] as a convex optimization program, as follows. 
First, let us rewrite (JHjl using the change of variables 7^ = 

T* — min YZ* fi 

{7i} 

s.t. (D-el)r-AghO 

P^^li^pT 1 , i = l,...,n, (7) 

where T — diag (7, ) . Therefore, the feasible set is convex in 
the space of variables 7,, i = 1, . . . , n. Furthermore, we now 
verify that the cost function Yl7=i fi (% rl ) ^ s a ^ so convex 
under Assumption 1 by computing its second derivative, 

^ E h (7-) = fi' (7- 1 ) \ t + 2/1(7-) ^> 0, 

where the last inequality is obtained from Assumption 1, 
taking into account that 7 4 _1 = Pi. 

The convex optimization program in |7]i allows us to 
efficiently find the cost-optimal allocation of vaccines to 
control the spread of an epidemic outbreak in a given 
contact network. In the following subsection, we illustrate 
our approach in a real social network. 



D. Numerical Results 

We illustrate our results by designing the optimal distri- 
bution of vaccines in an online social network when the 
cost vaccination function follows Q. We consider a social 
network with 247 nodes, and assume that the individuals in 
the network present the same recovery rate, <5; = 5 = 0.1. 
In this case, we can rewrite Q as a convex program 
with a convenient structure, as follows. First, defining a = 

\JT l -fil 1 ) , we have that 

E ^ (A) - a pi 1 - « E ft 1 = « Trace ( r ) - & > 

i i i 

where b = c*X^A rl - Hence, minimizing Ylifi(Pi) ^ s 
equivalent to minimizing Trace (r). Thus, the optimization 
problem in (|7]i can be written as the following semidefinite 
program (SDP): 

T* = min trace (r) 

r 

s.t. (S-e)T-AghO 

Pi 1 < 7i < §j\ i = l,...,n, (8) 

Given our network with 247 nodes, we now compute the 
optimal distribution of vaccinations in several cases. 

The network under consideration has a maximum eigen- 
value Ai (Ag) — 13.52. In our simulations, we choose the 
value of Pi to induce instability of the disease-free equilib- 
rium in the absence of vaccination. According to Proposition 
[T] if we had a constant infection rate /3j = /3 satisfying 
ft > ft c = 5/X\ (Ag) = 7.4e — 3, the disease-free equilib- 
rium is unstable. Hence, we choose a natural infection rate 
fti — ft > ft c to induce an unstable infection in the absence 
of vaccinations. In our simulations, individuals have the 
same natural infection rates fti = ft, and study three cases: 
ft € {1.2/3 C , 1.8ft c , 2.4/3J. We choose the value of ft. < ft c 
to induce a stable disease-free equilibrium in the case of full- 
force vaccination, i.e., we saturate all the individuals with 
vaccines to shift their infection rates to ft.. In our simulations 

— ^ _ 

we use a minimum infection rate ft. = 0.2/3^ = 0.2/3; 
hence, we obtain that ft. € {0.24^,0.36^,0.48^}. In 
other words, our vaccine reduces the infection rate to a 20% 
of the natural infection rate. Using these parameter values, 
we run three simulations, each one with a different ft. 

The results of our simulations are summarized in Fig. 
[2] Each one of the subplots in this figure correspond to 
a different value of ft € {l.2ft c , l.8ft c , 2.4/3J. For each 
value of ft we present a scatter plot with 247 data points 
(as many as individuals in the network), where each point 
has an abscissa equal to fi (fti) (the cost of vaccinating 
node i with optimal fraction fti) and an ordinate of di (the 
degree of i e V(G))- We observe that there is a strong 
dependence between the cost of vaccinating a node and 
its degree. In particular, we observe that there is almost 
an affine relationship between the vaccination cost and the 
degree of a node, with a saturation at the extreme cost 
values, and 1. Also, we observe that, as we increase the 
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Fig. 2. Vaccination costs versus degree in a social network with 247 nodes. 



value of j3 the vaccination costs tend to increase. This is Problem [2] can be stated as the following optimization 
because for larger ft the more virus is more infectious. problem: 



IV. Combinatorial Resource Allocation 

In this Section, we consider a combinatorial versions of 
the fractional vaccination problem studied in the previous 
section. In the fractional vaccination problem, the optimal 
distribution of vaccines is allowed to be in the feasible 
interval ft £ [ft, ft]- In the combinatorial vaccination 
problem, we restrict the resources to be in the discrete 
set, ft £ {ft, ft}. For this case, we propose a greedy 
approach that provides an approximation to the optimal 
combinatorial solution. We also provide quality guarantees 
for this approximation algorithm in Subsection |IV-B| The 
combinatorial vaccination problem can be stated as follows: 

Problem 2: Given a curing rate profile, {Si : i £ V (Q)\, 
and a vaccination cost function fi (ft) for ft £ |/3.,ftl, 
find the optimal distribution of vaccines to control tne 
propagation of an epidemic outbreak with an asymptotic 
exponential decaying rate e at a total minimum cost. 
The optimal distribution of vaccines in Problem [2] can be 
characterized by the set of individuals Ic C V (Q) that are 
chosen to be fully immunized, i.e., the infection rates are 
switched from ft to ft < ft for i £ Ic- Let us assume that 
the vaccination cost function takes the values fi (ft) = 
and fi (fij) = c i- These extreme values are achieved using 
the following affine cost function 

p i a \ A A ~ ft 
Ji (Pi) — c i n _ 7, ■ 
P_i ~ Pi 

Hence, the total cost of vaccination satisfies 

n 

^2 fi (A) = a c ^2 Ctfti - b c , 

i—l i 

where we have defined the constants ac — — ft) 1 
and be = ac 2~2i c iPi- Thus, since ac < 0, the optimal 
allocation of vaccines that minimizes 2~27=i fi (ft) ' s ^ e 
same as the one that maximizes ^ Cjft. Therefore, defining 
the vectors c = (ci, . . . , c„) T and b = (ft, ... , /3 n ) T , 



T% = max c T b 
C {ft} 

s.t. Ai (BA g — D) < —e (9) 
ft £ |^ 4 ,ft} , i = 1, ... 

The solution to this problem is combinatorial in nature. In 
the following subsections we provide a greedy approach that 
approximates the combinatorial solution, as well as a quality 
guarantee of our approach. 

A. Greedy approach 

In this subsection, we provide a greedy algorithm that 
iteratively updates the set of nodes that will be (fully) 
vaccinated in order to control the spreading of an epidemic 
outbreak. In each step of our algorithm, we denote the 
set of nodes that are chosen to be part of the vaccination 
group at St- We iteratively add to this group the node 
that provides the most benefit per unit cost, where the 
benefit of vaccinating a is the increment it induces in 
Ai (BAg — D). More formally, given a vaccination group 
St, we define the diagonal matrix of associated infection 
rates as Bs t — diag (ft) — (ft — /3 .)diag(ls t ), where ls t 
is the n-dimensional indicator vector for the set St. Thus, 
the benefit per unit cost of adding node i to St is measured 
by the function 

A( . g t) a ^ ( B * A s - D ) - A i ( B s t +{j} A e - D ) 

A conventional greedy approach could be defined by the 
iteration ft +1 = S t + {it} with ft = {} and i t = 
argmaxj A (i, Si), where this iteration is repeated until 
Ai (Bs t Ag — D) < — e is satisfied. Notice that the result- 
ing vaccination group is feasible and satisfies the spectral 
condition needed to control the spreading of an epidemic 
outbreak. 

In practice, we observe that a modification of this greedy 
approach provides better results. In this modified version, 
we start with a vaccination set ft = V (Q) (i.e., all the 



Parameters 


Metric 


Greedy 


Reverse Greedy 


Degree Threshold 


Centrality Threshold 


D* 


P = 2.4/3 c 


c'b 


3.6298 


3.6440 


3.2892 


2.4518 


3.9425 


P = 0.2/3 


XiiSB- 1 - A) 


0.0054 


0.0355 


0.0422 


0.1982 


n/a 


p = 1.8/3 C 


c'b 


3.0098 


3.0098 


2.9246 


2.0092 


3.1406 


P = 0.2/3 


Xi(5B~ l — A) 


0.0850 


0.1383 


0.0774 


0.2575 


n/a 


P = 1.20c 


c'b 


2.1484 


2.1484 


2.1201 


1.7369 


2.1787 


/3 = 0.2/? 


X^SB-'-A) 


0.4383 


0.4383 


0.6278 


1.0101 


n/a 



Fig. 3. Table with values of the objective function c'b and the residual value of Ai (<5_B 1 — Ag) for each possible value of ft. 



individuals are vaccinated) and iteratively remove individ- 
uals according to the iteration S t +i — S t — {jt} with 
jt = argminj A (J, St\ {j}), where this iteration is repeated 
until Ai (Bs t Ag — D) > — e is satisfied. The final vaccina- 
tion group is chosen to be St-i- Notice that, the resulting 
vaccination group is feasible and Ai [Bg t _ 1 Ag — D) < —s. 
We denote this approach the reverse greedy algorithm. 

Since our approach is heuristic for a combinatorial prob- 
lem, we provide a quality guarantee via Lagrange duality 
theory in the following subsection. 

B. Quality Guarantee 

Using Lagrange duality theory, we provide quality guar- 
antees for the performance of our greedy approach by 
computing the dual optimal D* c . 

Theorem 4.1: Given the optimization problem 

T c = max,, c T b (10) 
s.t. (D - el) B^ 1 - Aq>0 

A e {£.,&}, Vt, 

the primal optimal T c can be upper bounded by D* c 
computed according to the Lagrange dual 

min^.u l T u — trace(AgZ) (11) 



Dt 



s.t. Ui > CiPi + -^Z H Vi 

Pi 

Si 

Ui > Cij3. + — Zy Vi 

which is a convex Semidefinite Program. 

Proof: Notice that matrix in the semidefinite constrain 
can be written as (D - el) B^ 1 - Ag = J2i e i e 'i ~ Ag, 
where is the unit vector in the standard basis. From ( fT0| i, 
we construct the Lagrangian 



C{b, Z) = c T b + trace I Z I V e % e[4- - V: 



(12) 



where /3j € {/3., fti} is kept as a domain constraint and Z y 
0. See Section 5.9 of [20| for further details on the Lagrange 
dual of semidefinite constraints. Using the properties of trace 
to simplify and decouple we get 



£{b, Z) = ' s ^ j iciPi + -7fZ u \ - trace(ZA, 



(13) 



The dual objective is derived by maximizing the Lagrangian 
with respect to the primal variables 

Si . 



q(Z) = ^ ( maxci/3i 



Pi 



trace(ZAg). (14) 



Due to the decoupling in ( fT3| > the primal optimization in 
( fl4"] ) can be done for each node, independently. Since each 
node has only 2 options we can consider each case explicitly 
by defining 



max < dPi + yrZu, dP 



s A 

Pi 



(15) 



It is possible to compute m as a threshold function of Zu, 
but for the purpose of constructing the dual it is better to 
use an epigraph formulation to rewrite ( fl4| as 

q(Z,u) = J]] — trace(ZAg) (16) 



with the addition constraints that 



Ui 



> 



> 



CiPi + TrZu 
Pi 

c iP . + —Za. 

Pi 



(17) 



(18) 



Since the dual is a minimization and q(Z, u) is strictly 
increasing in u, either ( fl7| i or ( fT8] l must be achieved with 
equality, ensuring that the definition ( fl5| l is satisfied at 
the optimal point. To conclude, our dual (JTTJ is given by 
minimizing ( |X6] > subject to the domain constraint Z >z and 
the epigraph constraints ( fTT] ) and ( fT8j ). This is a standard 
form SDP as defined in section 4.6 of E0l . The solution 
Dq is guaranteed to satisfy D* c > Tq by weak duality, 
Il20l Section 5.2. ■ 
Theore m |4.1| tells us that for any optimization problem of 
the form (|10| we can get an accuracy certificate 



T* 



c T b <D%- c T b 



(19) 



by solving the dual ( fTT) . Since we do not have a strong 
duality, we do not expect c T b = D* c to be attainable (i.e, 
P*c < D* c ). 

Remark 4.1: The solution to the dual gives us some 
insight into the primal optimizers via the threshold solution 
to (BJ, 



Ui(Zu) 



CiPi ~t~ p Za 



if Zu < 



iMi 

if Z tl > fJip. 



(20) 



Relationship between Vaccination Strategy and Degree 




node degree, d 

Fig. 4. the relationship between the outcome of each algorithm and the 
degrees of the nodes. We represent degrees in log scale versus the fraction 
of nodes of a particular degree that are vaccinated in the solution generated 
by each algorithm. 



Relationship between Vaccination Strategy and Eigenvector Centrality 
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Fig. 5. Relationship between the outcome of each algorithm and the 
eigenvector centrality. In the abscissae, we represent the cumulative fraction 
of nodes with centrality greater or equal to a given value being vaccinated 
in the outcome of each one of the four algorithms under consideration. 



It appears we can deduce the primal optimizers b* from 
Z*, but in practice for most nodes i, Z* { = f3i/3.Ci/5i 
making it impossible to determine ft*. In some cases there 
are nodes that have not equal to the threshold. These 
nodes have their optimal action specified by Z*j and ( |20] >. 
This at least allows for a reduction of the dimension of the 
primal problem which due to its combinatorial form could 
be a very large improvement. 

C. Numerical Results 

Several papers in the literature advocate for vaccination 
strategies based on popular centrality measures, such as the 
degree or eigenvector centrality ETI . In this subsection, we 
compare our greedy heuristic to vaccination strategies based 
on centrality measures. In our simulation, we use the adja- 
cency matrix with 247 nodes previously used in Subsection 
III-D| and the same values for the parameters <5j = 5 = 0.1, 
ft - 6/X max _(A g ) = 7.4e - 3, ft e {1.2ft, 1.8ft, 2.4ft} 
and/3 = 0.2ft for all i. In Table pi we include the values 
of the objective function c'b ana the residual value of 
Ai (SB^ 1 — Agj for each possible value of ft. In each 
case, we run the greedy algorithm and the reverse greedy 



algorithm (both proposed in Section IV-Ai, as well as two 



previously proposed algorithms based on the degree and the 
eigenvalue centrality metrics. In the last column of Table [3 
we also include the upper bound provided by Theorem |44 
Observe that our greedy algorithms are always within 10% 
of the upper bound D* c . Furthermore, the reverse greedy 
algorithm is outperforms the others, specially those based 
on centrality measures. 

In Fig. |4j we illustrate the relationship between the 
outcome of each algorithm and the degrees of the nodes. 
In the abscissae, we represent degrees in log scale, and in 
the ordinate we provide the fraction of nodes of a particular 
degree that are vaccinated in the solution generated by each 



algorithm. We observe how all four algorithms completely 
vaccinate the set nodes with degrees beyond a threshold. 
On the other hand, in the range of intermediate degrees, 
we observe that degree alone is not sufficient information 
to decide the vaccination level of a node. In other words, 
simply vaccinating nodes based on degree does not always 
provide the best results. 

In Fig. [5] we illustrate the relationship between the out- 
come of each algorithm and the eigenvector centrality. In the 
abscissae, we represent the cumulative fraction of nodes with 
centrality greater or equal to a given value being vaccinated. 
We observe how all four algorithms completely vaccinate 
the set nodes with the highest centralities. However, since 
the curves in this figure are not monotonically increasing, 
there must be cases in which lower centrality nodes are 
vaccinated, but other nodes with higher centrality are left 
unvaccinated. In other words, vaccinating higher centrality 
nodes does not always provide the best results. 

The reason neither degree nor centrality adequately cap- 
ture the importance of nodes is that the eigenvectors of 
the matrix SB^ 1 — Ag change when the set of vaccinated 
nodes change. The shift in the eigenvectors is a result of 
the fact that optimal vaccination strategy actually depends 
on the parameters (3 and ft not just the network. With this 
in mind, we cannot expect an optimal solution to arise from 
an algorithm that depends only on the graph structure. Our 
algorithms work because they are greedy with respect to this 
effect. 

V. Conclusions 

We have studied the problem of controlling the dynamic 
of the SIS epidemic model in an arbitrary contact network by 
distributing vaccination resources throughout the network. 
Since the spread of an epidemic outbreak is closely related 
to the eigenvalues of a matrix that depends on the network 



structure and the parameters of the model, we can formulate 
our control problem as a spectral optimization problem in 
terms of semidefinite constraints. In the fractional vaccina- 
tion case, where intermediate level of vaccination are al- 
lowed, we have proposed a convex optimization framework 
to efficiently find the optimal allocation of vaccines when the 
function representing the vaccination cost satisfies certain 
convexity assumptions. In the combinatorial vaccination 
problem, where individuals are not allowed to be partially 
vaccinated, we propose a greedy approach with quality 
guarantees based on Lagrangian duality. We illustrate our 
results with numerical simulations in a real online social 
network. 
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